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ABSTRACT 

Dark matter substructure around nearby galaxies provides an interesting opportunity 
for confusion- free indirect detection of dark matter. We calculate the boost over a 
smooth background distribution of dark matter for gamma-ray emission from dark 
matter self-annihilations in tidal structure in M31, assuming a cross-section inversely 
proportional to the relative velocities of the dark matter particles as proposed by the 
Sommerfeld effect. The low velocity of the material in the structure results in a signif- 
icant increase in gamma-ray emission compared to both the background halo and the 
predicted emission for a velocity-independent cross section. We also calculate the ex- 
pected signal for Ferm% for reasonable choices of the dark matter parameters. We find 
that for a cross section proportional to the enhancement to the annihilation rate 
is sufficient to test the velocity dependence of the cross section by spatial correlation 
with the stellar component of the stream, given sufficient detector sensitivity. 

Key words: astroparticle physics-galaxies:individual:M31-galaxies:kinematics and 
dynamics-cosmology: dark matter-gamma-rays: galaxies 



their edges (Figure [T] left panel and right top panel). If the 
three features were all indeed produced in the same minor 
merger, the extremely high eccentricity required for the orbit 
of the progenitor implies that the 'shelves' are in fact radial 
fold catastrophes, otherwise known as caustics. This theory 
explains the sharp edges of the shells as the point where 
in-falling and outgoing streams of material stripped from 
the progenitor pass each other near the outer radial turn- 
ing point of their orbits. Since the motion is nearly radial, 
the projection of phase space into the r — Vr plane (Figure 
[l] right bottom panel) contains nearly all the information 
about the dynamics of material in the stream and shells. A 
caustic occurs at each point where the phase space stream 
becomes vertical in this projection, and the various features 
can thus be placed in chronological order of formation. This 
theory also predicts that near the edge of each shelf, the 
density will be significantly enhanced, as particles 'pile up' 
near the radial turning points of their orbits. 

Inferring the phase space distribution of the material 
in the shells and stream from this N-body model also re- 
veals that the relative velocity of material in the features is 
extremely low, especially in the tidal stream and the very 
edges of the two caustics (Figure [2|. At the caustic surface 
and in the stream, the local relative velocities can be less 



1 INTRODUCTION 

Like our own galaxy, the nearby Andromeda galaxy (M31) 
exhibits a wealth of stellar features with an accretion ori- 
gin, including dwarf galaxies, tidal streams, and a complex 
outer halo structure (j McConnachie et al.||2009 ). Unlike our 



galaxy, however, we view M31 from the outside, which in 
some cases facilitates the study of these substructures. Many 
of these tidal features are potential tracers of coUisionless 
dark matter associated with their progenitors, which under- 
goes similar dynamics to the stars. These features usually 
have distinctive, asymmetric shapes at large radii from the 
center of their host. This paper considers whether dark mat- 
ter in tidal structures in M31 could provide a possible indi- 
rect detection by ultra-sensitive gamma-ray observations in 
a confusion-free region around this nearby galaxy. 

One particularly prominent feature around Andromeda 
is a giant tidal stream that extends nearly radially away from 
the center of M31, commonly known as the 'Giant Stream.' 
This feature was first observed by Ibata et al. ('2001) and 



has since been studied in great detail. An N-body model 



of the stream by Fardal et al. (2006) has tentatively con- 



nected this stream with two other tidal features closer to 
M31's disc, known as the west and northeast 'shelves' be- 
cause of a relatively abrupt drop in surface brightness at 
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Figure 1. Two tidal features noted around the Andromeda galaxy (left panel, green and red points; top right panel, green and red 
dashed lines) correspond to fold catastrophes, or caustics, in individual phase- wraps of material from a dissolved satellite galaxy on a 
nearly radial orbit (right bottom panel). For this reason we refer to the shells in this work as Caustic 1 (green) and Caustic 2 (red), in 
the order in which they were formed. Another prominent tidal feature, the giant stream shown in orange, is the first structure to form 
in the merger, and hence is labeled '0'. 



Table 1. Mean and minimum estimated velocity dispersions in 
the features shown in Figure^ 



Feature (colour in Figure [l| km/s km/s 



Giant Stream (orange) 24 3.6 41842 

Caustic 1 (green) 70 7.3 29547 

Caustic 2 (red) 84 18 12263 



than 10 km/s (Table [T]). This is a result of the increasing 
thinness of the stream in phase space as time passes, an ef- 



fect sometimes known as 'gravitational cooling' (Mohayaee 
k Shandarin||2006t. 



Features of this type, though expected to be fairly com- 
mon, are as difficult to detect in the Milky Way as they are 
straightforward to find in sufficiently deep images of exter- 
nal galaxies. The sharp-edged shells seen in the star-count 
map of M31 would, when viewed from within Andromeda, 
look instead like large amorphous clouds spread over a huge 
fraction of the sky (Figure [3|. From this vantage point much 
more information about the phase space structure of the de- 
bris would be necessary to determine that the shells existed, 
whereas when viewed externally in a suitable projection the 
sharp edges immediately imply a nearly radial orbit for the 
progenitor. Thus, the existence of such a structure in M31 
represents a unique opportunity to study a system with well- 
constrained dynamics, thanks to its distinctive morphology 
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Figure 2. Projected phase-space plot of the tidal debris with a < 
200 km/s, colour-coded by the local velocity dispersion estimated 
from the N-body model (estimation method described in Section 
|2|. The cut in velocity dispersion excludes mainly material near 
the centre of the halo. The coldest material is found at the edges 
of the shells and in the tidal tail. For display purposes, a random 
selection of one-tenth of the particles are plotted here. 
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Figure 3. The N-body model of the tidal shells in M31, now 
viewed in Aitoff projection from a viewpoint at M31's centre. 
From this perspective the shells' sharp edges are virtually indis- 
tinguishable. 



as viewed from outside, at the closest range possible without 
full six-dimensional phase space information for stars in the 
shells. 

The existence of cold, high-density regions at large radii 
from M31's centre makes the tidal features an interest- 
ing candidate for indirect detection of dark matter in the 
Sommerfeld enhancement framework, where the interaction 
probability is boosted at low velocities. The Sommerfeld ef- 
fect was first introduced to boost the dark matter annihila- 
tion signal in order to account for the PAMELA observations 
of positrons and HESS and Fermi observations of an un- 
predicted high energy lepton component in the cosmic rays 



( Arkani-Ham ed et a H2009| ICirelh Cline|20Tol iLattanzi fc| 
Silk||2Q09| [March-Russell k West||2QQ9t . The annihilations 
of a TeV SUSY WIMP can be boosted by a factor of or- 
der 1000, as needed to account for the observed signal. This 
interpretation has been criticised on several grounds. Ex- 
cessive gamma rays (inverse Compton) and high energy an- 
tiprotons would be produced in the inner galaxy and exces- 
sive radio synchrotron emission in the outer galaxy if the lo- 



cal cold substructure persists at all galactic radii (Borriello, 



Cuoco & Miele 2009). The weakness in this critique is that 
the substructure is likely to be a strong function of galac- 
tic radius. Decrease in substructure concentration factor at 
large galactic radius and the effectiveness of tidal disruption 
at small galactic radii weaken these constraints significantly. 
Slatyer et al.|( [20TT| have systematically explored this effect, 
and conclude that our uncertainty about the radial depen- 
dence of the substructure contribution means that no strong 
constraints can be drawn from comparing signals at differ- 
ent Galactic radii. Any additional information about the ex- 
pected size of this contribution is therefore quite important 
for determining the viability of the Sommerfeld model. 

The strongest constraint on Sommerfeld boosting has 
come from considerations of delayed recombination of the 
universe following the last scatterings of the cosmic mi- 
crowave background radiation at z ~ 1000. The survival 
of the model appears marginal for WMAP constraints, and 



cantly weakened ( Slatyer et al.|20TT ). Additionally, this line 
of argument assumes that Sommerfeld boosting, quenched 
at the local value required to account for PAMELA, HESS, 
and Fermi observations, applies in a regime where the dark 
matter is much colder {/3 = v/c ^ 10~^^) than in the local 
halo substructure where /3 ^ 10~^. This is a huge extrap- 
olation that may not necessarily be justifiable in terms of 
general particle physics models. We believe it is important 
to test Sommerfeld boosting in a regime much closer to the 
local environment, for example in our galaxy, in ultra-faint 
dwarf galaxies, and in M31. 

Substructure in M31 can provide precisely such a test 
if we use old stars as dark matter tracers. The proposed 
test involves Fermi imaging of very similar substructures 
(at least in terms of velocity dispersion) to those invoked 
locally for the PAMELA excess. Without the Sommerfeld 
enhancement, the material is insufficiently dense to produce 



Planck will soon greatly improve these limits (Galli et al 



2009). However, if the contribution of substructures to the 



local signal dominates, then these constraints are signifi- 



detectable signal ( [Sanderson &: Bertschinger 2010), but 
with the enhancement signals can be boosted by a factor 
of up to 10^~^ at velocities comparable to those achieved 
in the tidal debris in this example. These boost factors are 
similar to those expected from local dwarf galaxies (Ro bert-[ 
[son Zentner|2009 ). With this in mind, we calculate in this 
work the boost and signal in the Fermi band from this tidal 
substructure assuming that a dark matter component of the 
unbound substructure follows the stellar component, as an 
example of the kind of result one might expect from tidal 
debris for this class of dark matter models. This particular 
example has then the additional advantage of a distinctive 
morphology that could allow it to be easily differentiated 
from a smooth dark matter halo. It also occupies an in- 
teresting niche between the bound substructures thought 
to dominate the signal in the Milky Way's outer halo and 
the more diffuse tidally disrupted substructure that [ Slatyer [ 
[et al .l propose contributes to the extragalactic gamma-ray 
background. In Section [2j we describe the method by which 
the phase space distribution in M31's halo and tidal sub- 
structure were modelled, and the results of tests for possible 
bias in our numerical methods. In Section [3] we present re- 
sults for the boost factor over the smooth halo as a result 
of the tidal substructure for different regimes of Sommerfeld 
boosting. In Section [4] we present maps of the flux in the 
Fermi band for two choices of dark matter model and show 
how those results may be scaled to other parameter choices. 
In Sections [5] and [6] we discuss the results and indicate paths 
for future work. 



2 MODELLING 

The rate F at which dark matter self-annihilations occur 
is proportional to the volume integral of the total squared 
number density of dark matter nJot , weighted by some func- 
tion S{v) of the relative velocity of particles whose form de- 
pends on the class of theories being considered. In our model, 
there are two distinct density distributions that contribute 
to the total density: the smooth halo distribution Uh and the 
tidal structure Us . The total rate can thus be separated into 
three different contributions for ease of calculation: one from 
interactions between dark matter particles in the smooth 
halo (denoted with a subscript hh), one from dark matter 
in the tidal structure interacting with dark matter in the 
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smooth halo (denoted with a subscript hs), and one from 
dark matter particles in the tidal structure interacting with 
each other (denoted with a subscript ss): 



Ttot 



J nl,Siv)dV 
J nlS{ah)dV + 2 

Thh + Ths + Fas 



/ 



nhnsS{vs)dV + 



/ 



nsS{as)dV 

(1) 



Here we have suppressed the position-dependence of the ar- 
guments for brevity, and denoted the volume element as dV. 
The argument to S{v) varies for these three terms. For Fhh, 
the correct relative velocity is the velocity dispersion an of 
the halo. For Fhs, the halo is assumed to have zero mean ve- 
locity relative to the debris, so the mean velocity Vs of shell 
particles is used. For Fss, the velocity dispersion of particles 
in the tidal debris, as is used. The shells and tail have such 
a low velocity dispersion that this last term is anticipated 
to dominate. 

To represent the phase space distribution of the ma- 
terial in the tidal shells, we used the N-body model con- 
structed by iFardal et al.l (|2006| 120071) to match the stel- 



lar component of the tidal debris. We assumed that the 
dark component tracks the stellar component and is of equal 
mass. These assumptions are admittedly an oversimplifica- 
tion but provide a good starting point for two reasons. First, 
the dark matter components of dwarf galaxies are thought 
to be more extended and less concentrated than the stel- 
lar component ( Penarrubia, Navarro McConnachie|2008a| ) 
but the starting conditions for the N-body model locate the 
satellite deep in the potential of M31, by which point this ex- 
tended dark halo would have been tidally stripped already, 
leaving only the dark matter within a tidal radius consistent 
with the stellar extent of the satellite. Within this radius, 
the dark matter is thought to contribute roughly equally 
with the stellar component to the potential of dwarf galaxies, 
providing some support for assuming comparable masses for 
the two components. This particular N-body model admits a 
dark matter component up to 2-3 times the total stellar mass 
(Fardal, private communication). Second, given the similar 
initial conditions of the stars and dark matter, the stellar 
shells can provide a visible starting point and template for 
searches because the radii of the shells, though perhaps not 
identical, will be similar. A justification for this is seen in 
simulations that follow both dark and stellar components of 
a nearly radial merger, as show n in Figure[4|(see Penarrubia, 



McConnachie & Navarro||2008| [Penarrubia, Navarro &: Mc- 



Connachie 2008b for full details of the simulations). The 



star shells, although formed well after the formation of the 
dark matter shells, clearly trace the dark matter shells and 
form at similar radii. 



2.1 Smooth halo background 

The halo used to model the tidal structure was also used to 
calculate both the smooth background signal and the signal 
from interaction between dark matter in the halo and pu- 
tative dark matter in the tidal shells. To make the N-body 



model, only the density profile p{r) was necessary (Geehan 
et al.|200"6 ); it is represented by the NFW relation (Navarro, 



Frenk &: White|1996 ) with scale radius rn and scale density 
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Figure 4. A recent simulation of a satellite falling into a Milky 
Way type galaxy (see [Penarrubia et al.||2008b[ for details of the 
simulation), showing the creation of shells of both stars (green) 
and dark matter (red) in a minor merger with a nearly radial 
initial orbit. The figure shows that star shells trace very closely 
the dark matter shells and despite the difference in initial distri- 
bution, the dark matter and star shells are at similar radii. The 
inset is a zoom on the satellite which is extremely resilient to the 
tides. 



ph,o, with the addition of a small core radius rcore to produce 
a finite central density: 



Ph(r) 



Ph,0 



(2) 



[(r + rcore) /rh][l + (r + rcore) /vh]^ '' 

as shown in the left panel of Figure |5] Note that this 
halo, used consistently for both the dark matter annihilation 
background and the N-body model, has a concentration of 
c = 25.5. This value is significantly higher than the concen- 
tration of a typical simulated, M31-sized dark halo grown 
from cosmological initial conditions, which in most cases is 
in the range 8-16 ( Navarro et al.||20i0| ). This high concen- 
tration is a byproduct of the dynamical model and relates 
to the uncertainty in the mass of M31's disc (the so-called 
'disc-halo degeneracy'); in this work it has the additional 
effect of producing a larger than usual background signal 
from the smooth halo density distribution in the innermost 
regions of the halo. 

The velocity dispersion a{r) of the halo was inferred by 
analogy with the high-resolution numerical studies of the 
phase space structure of cosmological haloes by jNavarro] 
et al. (2010). These studies confirmed the pseudo-phase- 
space-density scaling relation p/a^ oc r~^^/^, proposed by 



Bertschinger ( 1985) as a universal relation for cosmological 
dark matter haloes, over more than four orders of magni- 
tude in radius. The density profile of the halo was used to 
determine the radial profile of the velocity dispersion in the 
halo, scaled to a maximum velocity dispersion determined 
by rescaling one of the haloes studied in [Navarro et al.| The 
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Figure 5. Left: Density profile of the smooth dark matter halo used in this work. The grey solid lines highlight regions where the power 
law index is -3 (at large radius), -1 (at intermediate radius), and (at small radius). The long-dashed lines indicate rcore, the short-dashed 
lines indicate r/^, and the dot-dashed (green online) and dotted (red online) solid lines are the approximate radii of Caustics 1 and 2, 
respectively. Right: Velocity dispersion profile of the halo used in this work (solid black line). The profile peaks at rmax (long-dashed 
lines), which is just less than the scale radius (short-dashed line). The velocity dispersion of halo material in the region of the two 
caustics (noted as in left panel) is more than 10 times larger than the velocity dispersion of material in the caustics. 



average mean velocity of particles in the halo is assumed to 
be zero, and the resulting velocity dispersion profile is 



(Th{r) - an 



p{r) 



1/3 



5/8 



(3) 



where it can be shown that 



and 



729 



Pmax Y792^^ ^ 0-41p/x (4) 



The velocity dispersion profile is shown in the right panel of 
Figure [5] 

Using the analytic expressions for p and a, we then 
compute the quantity Fhh analytically: 



ml J 



dVplSiah) 



(5) 



with ph given by Equation ^ and an given by Equation 
([3|. For consistency with the A/'-body model, we normalise 
the expression to the number density of A/'-body particles by 
dividing by m^. The integral in Equation ([5| could be taken 
over the entire volume of the simulation to estimate the total 
halo flux, but this value depends very strongly on the choice 
of a core radius for the halo, which is not constrained by 
the dynamical model. This sensitivity, however, is confined 
to a tiny region right at the centre of the halo, equivalent to 
the central four pixels or so for Fermi. To avoid the strong 
dependence on a parameter that is so ill-constrained, we 
instead compare the signal from the halo to that from the 
tidal structure pixel- by-pixel over the field. Equation ([s} is 
evaluated for each pixel separately over the entire line of 
sight (z) and an area on the sky (Ax, A^) corresponding 
to the resolution of Fermi (about 0.1 degrees, or about 1.4 
kpc at the distance of M31). The boost is calculated by 
comparing the signal in each pixel from the halo and tidal 
structure. The uncertainty about the core radius leads to 
unreliable estimates of the boost only for the central few 
pixels, while the region of interest is at larger radii where 
the mass profile is somewhat better constrained, and where 
the boost factor does not depend so strongly on the slope of 
the mass profile. 



2.2 Tidal structure 

The tidal structure, including two shells and the giant 
stream, is represented by an N-body realisation based on 



the model constructed by [Fardal et al. (2007). Although 
this model is by no means the single best fit to the available 
data, it is at least a local best fit that provides a plausi- 
ble dynamical origin for the debris and a valuable tool for 
inferring the phase- space structure. Numerical methods are 
necessary to estimate the signal from the tidal structure; 
here we describe these methods and their limitations. 

The integrated squared density, weighted by a factor of 
either 1/v or was estimated using the optimal proce- 

dure identified in |Sanderson &: Bertsching"er| ( |2010| ), with the 
addition of estimates for the mean velocity (used to account 
for interactions between shell and halo dark matter parti- 
cles) and the velocity dispersion for material in the shell. 

Moments of the velocity are calculated as follows. First 
the mean velocity v at the centre of the current Riemann 
volume, located at position x, is estimated using 



v(x) 



(6) 



Here and in the following, the hat symbol indicates an es- 
timator that recovers a smoothed field from the discrete 
N-body representation. The average relative speed Vs(x) is 
then calculated by taking the magnitude of the mean veloc- 
ity vector: 



'Os(x) — Vv - v. 



(7) 



The quantity Vs(x) represents the relative velocity between 
material in the shell and material in the halo at point x. 

The mean velocity at the point x is used to compute 
the nine-component, symmetric velocity dispersion tensor 
cr*-^, for the orthogonal directions {i,j} G {x,y,z}, at the 
same position x: 

n=l 

The average one-dimensional velocity dispersion crj(x) 
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is calculated by summing the three eigenvalues of the 
velocity dispersion tensor estimated with a^^ (x) (the lengths 
of the orthogonal axes of the velocity ellipsoid) : 



(9) 



We determine as by simply taking the square root. This 
quantity represents the relative velocity between particles 
in the shell at point x. 

Finally, the two components of F involving the A/'-body 
model are computed using the estimators 



f hs = — Vihs,iPh{^i)S{vs) 



Viev 



and 



(10) 



(11) 



where the Riemann sum is over the volumes Vi making 
up the target volume V, and hi and n'^i are estimated at 
the centre x^ of each Riemann volume as in [Sanderson fc| 
Bertschinger (2010). The halo density ph is evaluated at the 
centre of each Riemann volume for consistency with the po- 
sitional accuracy of the density estimates. Because the Rie- 
mann volumes are generally small compared to the gradient 
of the halo density profile in the regions of interest, the differ- 
ence between this method of evaluating ph and the analytic 
integral over each pixel used to calculate Fhh should likewise 
be small. 

Recovery of smooth fields from a discrete representation 
can be sensitive to various discreteness effects, including the 
choice of smoothing number Ns and resolution Np and the 
local gradient of the density, especially the existence of sharp 
edges in the distribution. Because of the complexity of the 
method for determining the velocity dispersion, we used nu- 
merical experiments to calculate the bias, variance, and rms 
error of as for several different values of Ns, and Np over a 
range of velocity dispersions. We looked for variations with 
these parameters, as well as those due to edge effects or 
density gradient (which is high near the caustic). 

For the purposes of the numerical experiments, we de- 
fine the expectation value of the estimator, E(as), as the 
mean of the estimated values of as over a given region of 
the sample: 

E{as) = -^XiVsuba,(x,), (12) 

t — 1 

where A^sub is the number of particles falling in that region 
of the realization and x^ is the position of the zth parti- 
cle, with as defined as described above. A "region" could 
be the entire sample, in which case A^sub = Np, but we also 
compared subsets that included and excluded edges or caus- 
tics. We compared E{ds) to the input value of the velocity 
dispersion, (Tin, by computing the bias B and variance V of 
E(as), defined as 



B = 



E{as) 



and V = ^ ^2 [o-s(x^) - E{as 



(13) 

With these definitions, B represents the average fractional 
systematic difference between the input and estimated val- 
ues of (Js, and similarly V measures the average fractional 



variation of the estimates from their mean (the "noisiness" 
of the estimator). The square root of V, often referred to as 
the standard deviation, is a measure of the spread of all the 
individual estimates of as. 

One can further quantify the performance of the estima- 
tor by combining B and V in the root-mean-squared (rmse) 
error, defined as 



rmse)^ = + V : 



Nsuhcrf^ ^ 
Hi I— I 



[(3■s(x^) - an 



(14) 



which includes both error from a noisy estimator (in V) and 
error from a biased one (in B). It is important to consider 
the relative contributions of B and V to the rmse, however, 
so we will discuss all three quantities below. 

We calculated the bias, variance, and rmse of the expec- 
tation values from random realizations of three-dimensional 
distributions with and without caustics. The realizations 
were generated at a range of resolutions between Np = 10^ 
and Np = 10^, with the number of particles in a given re- 
alization chosen from a Poisson distribution with mean Np. 
We used a range of velocity dispersions between 10~^ and 
10~^ and a time unit of 1 to create the caustics. The caustic 
width is simply 6x = at in our map, so this creates caus- 
tics with a thickness of 10~^ — 10~^ relative to the units 
measuring box size as shown in Figure [6] The same veloc- 
ity dispersion used in the caustic mapping was also used 
to assign random velocities to each particle: the velocities 
were reassigned to caustic particles for consistency with the 
uniform case, and in order to isolate the effect of density 
gradients on estimates of the velocity dispersion. The veloc- 
ity dispersion was calculated at the location of each particle 
in each sample at Ns = 10, 20, and 30. 

We tested for edge effects by comparing the estimates of 
as for two types of test distributions with a uniform veloc- 
ity dispersion: one with a caustic in it and one with uniform 
density. The estimator described above was used to estimate 
as at the location of each particle in a given sample. To test 
for possible effects of edges and density gradients, we defined 
two regions in each sample: one near the edge of the distri- 
bution and one in the centre of the sample (Figure [gJ. The 
edge region in the caustic distribution is aligned with the 
caustic surface to probe possible bias from the high density 
gradient in the caustic, and spans several times the scale 
width of the caustic. The sizes of the regions are adjusted 
so they all contain about the same number of particles. The 
velocity dispersion was estimated at the locations of parti- 
cles in a given region using all particles in the realization. 
If the edges or density gradient affect the estimation of the 
velocity dispersion, we expect to see a difference in the bias 
and/or variance of the estimates for regions near the edges 
or near the caustic compared to regions that exclude the 
edges. 

To illustrate the effect of edges, we present results for 
Ns = 10 in a pair of realizations with Np = 10^ (correspond- 
ing to A^sub 10^"^), since these values of Ns, Np, and A^sub 
are appropriate for the N-body realization of M31. We found 
that the rms error on estimates produced by as was domi- 
nated by the large variance resulting from the low value of 
A^s and was about 14 per cent in all cases. The bias, which 
would indicate a systematic error in the estimator, was also 
independent of the density gradient or the presence of edges 
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Figure 6. Regions for testing possible edge and density-gradient effects in the estimation of the velocity dispersion ds. The test compared 
different regions (edge and center) in two test distributions: one with a caustic along the x direction (left panel) and one with uniform 
density (right panel). The edge region in the caustic distribution spanned the caustic surface (red box) while the center region avoided 
all the edges and had a small density gradient (orange box). The edge region in the uniform distribution covered all the edges (half are 
shown in purple) while the central region avoided all of them (green). 
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Figure 7. The bias (filled points) and rms error (error bars are 
1/10 the rmse) of E(as) do not appreciably differ between the 
various regions of the sample, or between the samples as a whole. 
Print version: squares denote the sample with the caustic; circles 
denote the uniform-density sample. Online: The colours in the 
figure correspond to the regions depicted in Figure [6] magenta 
indicates the entire sample with the caustic and blue indicates 
the entire uniform-density sample. 



in the distribution (Figure [7|. This rms error corresponds 
to an uncertainty of about 25 per cent in the Sommerfeld 
coefficient, which is comparable to our uncertainty about 
the details of the phase space structure in this tidal debris 
and certainly less than the uncertainty about the particle 
physics model of the dark matter. 

The N-body model of the debris uses about 1.3x10^ 
particles to represent the entire tidal structure, of which 
about 4 X 10^ end up in Caustic 1 and 2 x 10^ end up in 
Caustic 2. This level of resolution is sufficient to resolve the 
density distribution of the material if a suitable estimator 
with a suitably small smoothing number {Ns = 10) is used 
(Sanderson & Bertschinger 2010), but this does not guaran- 
tee that the velocity structure of the material is adequately 
resolved. To maximise the resolution of the velocity struc- 
ture we would like to use the smallest possible smoothing 
number to estimate moments of the velocity distribution as 
well. 



In order to understand how the choice of smoothing 
number and the resolution of the N-body representation af- 
fected the sensitivity of the calculation to small velocity dis- 
persions, we computed the bias and variance for estimates of 
a in uniform and caustic density distributions with different 
input values of the dispersion atrue, at different smoothing 
numbers Ng, and at varying resolutions Np. We find that 
the estimator can reliably estimate velocity dispersions as 
small as lO"'^, with no indication that the bias is dependent 
on a (Figure [s] left panel). If the box size is rescaled to the 
approximate size of the caustics in the N-body simulation 
of M31, (Jin = 10~^ corresponds to caustics of width ^ 5 pc, 
and velocity dispersions of about 5 km/s. The caustics in 
M31 have widths closer to 1 kpc, so their density gradient 
is always many times smaller than those tested though the 
minimum velocity dispersions are comparable to this limit 
(Table [T]). Choosing a larger Ns does slightly reduce both 
the bias and the rms error, but will make the estimator less 
sensitive to small-scale changes in the velocity dispersion. 
Increasing the resolution also has only a small effect on the 
bias and virtually none on the rms error (Figure [s] right 
panel), which is again dominated by the variance. Based on 
these tests we conclude that the current level of resolution 
of the N-body realisation and the choice of A^s = 10 will 
recover adequately unbiased estimates of the velocity dis- 
persion, sufficient for the required level of accuracy in this 
work. 

3 BOOST FACTOR 

The boost factor is defined as the enhancement over the 
smooth halo provided by the tidal structure: 

^ ^ Fhs + Tss ^^^^ 
Thh 

The boost factor is independent of the normalisation of 
S, the branching ratio, and other quantities that are deter- 
mined by the particular particle physics model of the inter- 
action. Thus, it usefully isolates the effect of the different 
dependences of the cross-section on velocity without intro- 
ducing all the complexity of the parameter space of dark 
matter models. We considered two different power laws for 
S{v) motivated by previous studies of the Sommerfeld effect 
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Figure 8. The ability of the estimator to recover velocity dispersions is fairly insensitive to most numerical parameters. The average 
bias using Ng = (10, 20, 30) is shown in (darkest to lightest gray; blue, green, red in online version) in the left panel, for various input 
values of crtrue in kpc Myr"-*^. Increasing the smoothing number Ng (left panel) or the resolution Np (right panel) slightly improves the 
bias and rms error regardless of whether the density distribution includes a caustic (solid lines) or is uniform (dashed lines). In the right 
panel, Ns = 10 and Np is varied; above log^^Q Np =4.5 there is not an appreciable difference in the bias. In both panels, the error bars 
represent 1/lOth of the rms error of the estimator. 



e.^. Arkani-Hamed et al.|2009|[Bovy|2009|rLattanzi &: Silk 



2009||MardQn et al.|2009t :\s(i;) oc 1/v and S{v) ocl/v'^. For 
comparison we also consider the velocity-independent case 
S=l. 

Both velocity-dependent cases (Figure [9] center and 
right panels) provide a significant, position-dependent en- 
hancement of the tidal structure relative to the background 
and relative to the velocity-independent case (Figure left 
panel). As expected from Figure |2] the most significant 
enhancement is at the edges of the two shells and in the 
stream, where the density is highest and the velocity disper- 
sion is lowest. The enhancement compared to the velocity- 



independent case is highly non-linear (Figure 10) because 
of the correlation between the density and velocity disper- 
sion, which is a product of the phase-space streaming and 
the radial symmetry of the system's dynamics. 

Surprisingly, the tidal stream produces an enhancement 
that rivals or exceeds that of the shells. In retrospect, exam- 
ination of the right panel of Figure ^ and the right panel of 
Figure [2] shows that the tidal stream is just as cold as (and 
perhaps colder than) the shells, its material falling radially 
inward in a narrowly collimated and fairly dense band. The 
tidal stream is fairly dense at radii even larger than the edges 
of the shells, where the halo is of almost negligible density, 
leading to an even larger boost factor. 

Of particular interest is the prediction that if the cross- 
section to dark matter self-annihilation is proportional to 
1/^^, the emission from this tidal structure should be as 
bright as the background halo at the edges of the shells and 
in the densest part of the tidal stream. This finding can 
be used to test the velocity-dependence of the dark matter 
cross section if the halo of M31 is detected in gamma-rays 
by Fermi. If the halo is detected. Figure [9] predicts that 
zeroth-order departures from a smooth emission distribu- 
tion should be observed if the cross-section depends on l/v"^ 



or higher order, and that these departures should be corre- 
lated with the spatial distribution of tidal structures around 
M31. Likewise, if the cross-section depends on 1/v, depar- 
tures should be observed at the ten percent level, although 
this may be beyond the range of current instruments. These 
predictions are independent of the specific model of the dark 
matter particle, and are based solely on the assumption of 
a form for the velocity-dependence of the cross section. If 
no such departures are observed, the class of models with 
velocity- dependent cross sections of that form can be ruled 
out. 



3.1 Spatial correlations 

The particular morphology of the tidal features can signif- 
icantly improve the chances of a detection for a low sig- 
nal rate, by correlating the stellar map with the gamma-ray 
map. Here we demonstrate a coarse version of this by divid- 
ing the map into several regions, three centred on a feature 
and one without significant tidal boosting, and calculating 



the boost in each of these regions separately (Figure 11). 
This coarse contrast method can also allow for slight de- 
viations between the dark matter and stellar distributions, 
though our calculation assumes perfect tracking. In practice, 
the boosts would be calculated by fitting a smooth, spherical 
halo profile to the radially averaged observed distribution of 
gamma rays, which does not show much deviation for any 
of the interaction models we considered, and comparing the 
observed and fitted signal in each pie-shaped region inspired 
by the arrangement of the tidal debris, which is highly asym- 
metric. Including a region assumed to have no boost gives a 
built-in measurement of the sensitivity of the comparison. 

The most prominent feature appears to vary based on 
the interaction model that is used. For S oc l/t'^, the tidal 
tail is the most prominent, deviating from spherical symme- 
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Figure 11. Comparing the boost in different regions can reveal 
deviations from spherical symmetry in the map, even if the signal 
from the substructure itself is not fully resolved, for some inter- 
action models. The inset shows the regions over which the sum 
was calculated. 



try at the 1-percent level. For S oc the edge-on shell 
appears at the 0.5-percent level but the tidal tail is indistin- 
guishable. This is because the cross-term Fhs is larger than 
the shell-shell interaction Fss for this case, so that the struc- 
ture at the smallest radius (i.e. largest halo density) is the 
brightest. Without a velocity-dependent cross section, none 
of the features is distinguishable from the background. 



4 GAMMA-RAY SIGNAL 

In order to determine whether the test described in the pre- 
vious section could be performed with Fermi, we estimated 
the gamma-ray signal from the halo and tidal substructure 
for the two forms of S{v) described above. In this section 
we discuss the calculation of the signal and its scaling with 
various parameters. 



4.1 Calculation of the signal 



We follow the notation of Fornengo, Fieri & Scopel (2004) 
as adapted by Sanderson & Bertschinger|( |2010| to calculate 
the gamma ray signal, dN^/dt, for Fermi. As in Fornengo 
et al. ( 2004 ) , we separate the dependence on the phase-space 



distribution of material from most of the details of the dark 
matter particle model: 



nn 



dE^ 



-A{E^)dE^dE^ 



(16) 



The first term, $ , depends only on the mass and velocity 
distribution of the dark matter and the velocity-dependence 
of the Sommerfeld effect. This term is independent of the 
energy E^ of the gamma rays produced in the interaction: 



4^1 dWi^)SH^)\ 



(17) 



The signal depends on the distance, to M31 and on the 
local relative velocity v of the dark matter as well as the lo- 
cal mass density p. Of course, the rate at which dark matter 
particles interact with each other really depends on the num- 
ber density, not the mass density, but since the dark matter 
mass is model-dependent there is a corresponding factor of 
1/m^ in the second term, d^^ /dE^, which depends on the 
particular model of dark matter being used. This term de- 
scribes the spectrum of gamma rays produced for a given 
dark matter model: 



d^^ _ {av)o dN^ 
dEj 2m| dE^ 



(18) 



The cross section {av)o denotes the value of the cross section 
without Sommerfeld enhancement. dN^/dE^ is the spec- 
trum of gamma-rays produced in a particular dark matter 
model. 

The total signal in a given detector, dN^/dt, is calcu- 
lated by integrating the spectrum of observed radiation over 
the energy range of the detector, weighted by the effective 
area for the detector ; the only other detector-dependent 
piece is the lower limit of the integral ^min, the threshold 
energy for detection of gamma rays, chosen to be 100 MeV 
for consistency with the Fermi sensitivity range. (The upper 
limit of the integral is simply the dark matter mass, as re- 
quired by energy conservation.) For the Fermi LAT, whose 
effective area is roughly energy-independent above 1 GeV 
dFermi LAT Collaboration||2009| ), we can calculate the flux 



independent of the effective area, since we are interested 
mainly in determining whether the structure is above the 
detection threshold. 

With this simplification, we find that is simply pro- 
portional to the total yield N^{Emin) above ^min, so that 



^''{En 



{crv)o 
2ml 



N^{En 



(19) 



4.1-1 Sommerfeld enhancement 

[Arkani-Hamed et al.| ( |2009t , [Lattanzi Silk| ( |2009| ), and 

others have shown that the Sommerfeld interaction can be 
easily modelled as a Yukawa force with coupling constant 
a and a mediating particle mass mi. In this analogy, the 
solution to the radial Schrodinger equation with a Yukawa 
potential exhibits two characteristic behaviours of the cross 
section enhancement S = (av) / {av)o. At very low relative 
speeds, for resonant values of the mass ratio m^/mi, S is 
proportional to the inverse square of the relative speed: 



(20) 



where 13 = v/c is the relative velocity of dark matter parti- 
cles. A 'low' relative speed is small compared to the charac- 
teristic velocity 



/3* 



ami 



(21) 



' pixel 



Previous work on the Sommerfeld effect has proposed values 
of around 10~^ for the coupling constant a, and a wide range 
of m^ and mi . We choose a = 1/30 for this work, noting that 
the enhancement only depends on the quantity (3* and may 
thereby be rescaled to any combination of ami/m^. The 
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first resonant peak occurs when mi/m^ ^ 0.2 (Lattanzi & 
Silk|[2009l ) 



The N-body reaUsation gives velocities in the units kpc 
Myr~^ and hence introducing appropriate units and values 
leads to the expression 



kpc Myr~^ 



(10-2) 



mi/m^ 

10-2 



(22) 



where S'o^' = 9.4. 

Away from resonance, the enhancement takes on the 
form: 



13 Vkpc Myr-i 



(10-2) ' 



(23) 



where 6*0^ = 9.6. This form is also valid at intermediate 
values of /3 between /3* and a. 



4.1.2 Phase-space distribution factor 



With an N-body realisation, the integral in Equation (17) 
is calculated in terms of the number density n for N-body 
particles of mass m^, so that: 



Cres 2 

Dq nip 



(10-2) 



mi/r 



J (i^xn^(x) 



^(x) 



47rd2 VlO-27 V 10"^ J J ' [kpc Myr-i 

(24)" 

for the resonant interaction and similarly for the non- 
resonant regime. The number density n is likewise in units 
of kpc~^ and the volume element d'^x in units of kpc*^ and 
hence the result of the numerical integration returned by the 
density estimator is 



kpc~^ Myr^ 
on resonance and 



/ 



kpc^ 



kpc Myr 



/ 



kpc^ 



kpc 



kpc Myr 



kpc 



i;(x) 



kpc Myr 



(25) 



(26) 



off resonance. 

Using the distance d = 785 kpc to M31 and the particle 
mass nip ^ 10^ Mq from the N-body simulation, and con- 
verting to standard units for the quantity we obtain the 
phase-space-dependent factor: 



785 kpc 

2 



104 M, 



(10-2) 



kpc ^ Myr^ 



mi/m^ 

10-2 



(27) 



where 

nant interaction and 



^C,nr ^ ^C,nr 



1.75 X 10"^^ GeV^ kpc cm-^ for the reso- 
(1F2) 



785 kpc 

2 



104 M 



kpc"'* Myr 



(28) 



where $0 = 1.79 x 10"^^ GeV^ kpc cm-^ for the non- 
resonant interaction. 



4.1.3 Particle physics factor 

The particle-physics factor in the flux can be written as 



3 X 10-26 cm3 s- 



- ) i V 

1 7 vio Gevy 



(29) 

whereto =4.63x10 ^ cm^ kpc ^ GeV ^ s \ The gamma 
ray yield is usually of order 1 per annihilation or less, 
depending on the dark matter model ( Mardon et al.||2009 ). 



4.1.4 Complete expression 



Combining Equation (27) or (28) with and (29) gives the 
master equation for including all the scalings: 



1 ( \ 

-< ' llOGiV J 

mi/mx\ f a \ f {(Jv)q 
UO-2J 



10- 



d 



3 X 10-26 cm3 s-i 

2 / „.^„2 



^104 Mo y V 785 kpc 
for the resonant process, and 



kpc~^ Myr^ 



,0) 



;.3 X 10 



AT -2 -1 ( A"^ 

^-^^ ^ (lOG^J 
VlO-2y V3 X 10-26 cm3 s-i ) 



(31) 



, 104 Mo y V 785 kpc J y kpc-^ Myr ^ 
for the non-resonant process. The comparable expression 



without the Sommerfeld boost is ((Sanderson & Bertschinger 



2010) 



9.6 X 10"^^ A^^ cm"^ s"^ 
""llOGevJ 1^ 



104 M, 



3 X 10-26 cm3 s-i 
2 / 1 \ -2 



d 



785 kpc 



^(^^ 



kpc 



(32) 



4.2 Results for the Fermi band 



Using Equations ( 30 ) , ( 31 ) and ( 32 ) , we produce maps of the 



total flux in gamma-rays in the Fermi band, including 
both the halo and the substructure, for several scenarios. 
Figure [12] compares the results for resonant, non-resonant, 
and non- Sommerfeld cases for two different choices of m^'. 
10 GeV and 1 TeV. The former is optimistic but realistic for 
models with no Sommerfeld boost; the latter is characteris- 
tic for models with a Sommerfeld boost. We include results 
at 10 GeV for Sommerfeld-like boosts for completeness, al- 
though a particle model for such an enhancement at low 
nix does not exist to our knowledge. However, we do note 



that Slatyer et al. (2011) point out that mediating particles 



with masses even lower than a few GeV (the lowest con- 
sidered here) cannot be ruled out by current measurements 
given the uncertainty about the distribution of substructure 
in the Galaxy, so this panel may yet be relevant. 

The case S oc 1/v is brighter than 5* = 1, but only the 
halo is visible in both cases: the tidal features are below the 
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Table 2. Total flux from the halo and tidal structure, compared 
to the recent Fermi detection of M31. 



Source 


^^,m^ = 10 GeV, <l>^, = 1 TeV, 

— 9 —1 —9 —1 

cm ^ s cm s 


$NS 

^SS 


1.6 X 10-9 
3.8 X 10-12 
4.8 X 10-14 


1.1 X 10-12 
2.7 X 10-1^ 
3.4 X 10-1'^ 


^nr 

ss 


5.0 X 10-8 
8.2 X 10-11 
5.8 X 10-12 


2.4 X 10-11 
4.0 X 10-14 
2.8 X 10-15 


*hh 

^res 
ns 
<i»res 

^ss 


2.9 X 10-'^ 
1.5 X 10-10 
2.8 X 10-10 


1.4 X 10-10 
7.2 X 10-14 
1.4 X 10-13 






9 X 10-9 



^Fermi LAT Collaboration] ( |2010t 

smooth emission by several orders of magnitude. In Table 
|2]we see that the cross-interaction signal $hs is larger than 
the signal $ss from interactions within the debris, so that 
overall the signal from the substructure scales only linearly 
with the substructure density at leading order. However, 
the structure may still be marginally detectable using the 
coarser test described in Section |3.1| given sufficient sensi- 
tivity to detect the halo at the appropriate radii. 

In the case where S oc l/t'^, the enhancement is non- 
linear enough in both p and a that although the center of 
the halo is still the brightest part of the structure, the tidal 
features stand out above the halo at their radii. Table [2] 
shows that in this case $ss > ^hs, indicating that for this 
case the leading-order signal really scales with the square of 
the substructure density. This structure is still at least an 
order of magnitude below Permits current sensitivity regime, 
but if such a sensitivity were achieved, a search for deviations 
from spherical symmetry in the gamma-ray emission would 
be able to test the velocity-dependence of the interaction 
cross-section. 

We also note that the values of $hh in the table depend 
primarily on the inner slope of the mass profile as discussed 
in Section|2] so the fact that they exceed the measured signal 
from M31 is merely a reflection of the uncertainty of this 
parameter. 



5 DISCUSSION 

Thanks to their low velocity dispersion and relatively high 
density, cold tidal streams and young caustics can provide 
a significant boost to the dark matter self-annihilation rate 
if the cross section is non-linearly dependent on the relative 
velocity, as in the Sommerfeld scenario. The particular mor- 
phology of tidal streams, their location far from the centres 
of galaxies, and the apparent tracking of the stellar and dark 
components also make these features an attractive place to 
search for an annihilation signal, as the correlation with the 
stellar shape makes it easier to differentiate such a signal 
from a smooth halo distribution. In fact, at distances typ- 
ical of tidal debris 30-50% of the halo's mass may be in 
streams dMaciejewski et al.|2010| [Wang et"aL]|2011t . 

In this work we computed the boost to a smooth 



background from tidal substructure observed in star-count 
maps of the Andromeda galaxy for two different velocity- 
dependent cross sections proposed for Sommerfeld-type an- 
nihilations between weakly interacting, massive dark matter 
particles (in which a light helper particle forms in an inter- 
mediate state between the dark matter annihilation and the 
production of standard- model particles). We used an exist- 
ing N-body model of the structure to estimate the density 
and relative velocity of the material in the substructure with 
suitably unbiased algorithms and a reasonable choice of nu- 
merical parameters. We found that in both cases, the emis- 
sion from the tidal structure could locally boost the annihila- 
tion emission by up to factors of 5. The case where the cross 
section (av) (x produced a boost large enough for the 
signal from the tidal structure to outshine the smooth halo 
at large radii, though the estimated signal is several orders 
of magnitude below the current Fermi sensitivity for reason- 
able choices of the dark matter parameters. However, if an 
instrument with the required sensitivity existed, a search for 
emission from the tidal substructure would constitute a test 
of the velocity-dependence of the dark matter cross section, 
since only in the l/v"^ case is the tidal structure visible. 



6 FUTURE WORK 

One intriguing result of this work is that tidal streams that 
are sufficiently massive and collimated (i.e., young and/or 
cold) can produce significant local boosts of the annihila- 
tion signal if the cross-section is velocity-dependent. Such 
streams could be a significant contributor to the gamma-ray 
luminosity of the outer regions of haloes in the Sommerfeld 
scenario, especially if cosmological simulations accurately 
predict the percentage of streams. State-of-the-art cosmo- 
logical simulations of individual galactic haloes can resolve 
this coarse-grained phase space structure, and could be used 
to estimate this contribution. 

Low-frequency radio observations could also be used to 
search for signals from dark matter in tidal substructures, 
via channels that produce high-energy electrons that then 
give rise to synchrotron radiation through interactions with 
the galactic magnetic field. A map of the polarisation must 
be correlated with the observed stellar stream, under the 
assumption that the dark matter and stars track each other, 
to search for such a signal (Zaroubi, private communication). 
In future work, we will consider whether this channel could 
produce a signal detectable with a low-frequency array such 
as LOFAR, either in M31 or in high-latitude streams in our 
own Galaxy. 
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